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Abstract 

In this work, we investigate the intrinsic limits of subthreshold slope in a dual gated 
bilayer graphene transistor using a coupled self-consistent Poisson-bandstructure solver. 
We benchmark the solver by matching the bias dependent bandgap results obtained from 
the solver against published experimental data. We show that the intrinsic bias dependence 
of the electronic structure and the self-consistent electrostatics limit the subthreshold slope 
obtained in such a transistor well above the Boltzmann limit of 60mV/decade at room 
temperature, but much below the results experimentally shown till date, indicating room 
for technological improvement of bilayer graphene. 
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The excellent transport properties of single and multi-layer graphene hold promise to build 
ultra-fast transistors with excellent on state characteristics pQ-[lj. However, the lack of signif- 
icant bandgap in such systems has been one of the major roadblocks to achieve low off state 
current and hence high on/off current ratio |2J-|4J. Recently, it has been found, both theoreti- 
cally and experimentally, that a bandgap can be opened up in a bilayer graphene (BLG) using 
an external bias A recent experiment shows that a bandgap as large as ~0.3eV can be 

created in a BLG depending on the external bias [6] . This relatively large bandgap is promising 
to obtain low off state current, hence an improved on to off current ratio and subthreshold slope 



The aim of this paper is to discuss the intrinsic limits of the subthreshold slope in a BLG 
metal-oxide-semiconductor field effect transistor (MOSFET). We assume a long channel tran- 
sistor with an ideal gate dielectric interface. To achieve this, we develop a self-consistent 
Poisson-bandstructure solver for a dual gated BLG device. We have simulated the magnitude 
of the bias dependent bandgap tuning in a BLG system which matches closely with the pub- 
lished experimental data [6] validating the solver. We also show that the gate bias dependent 
electronic structure and the self-consistent electrostatics play a central role in determining the 
intrinsic limits of the subthreshold slope in a BLG transistor, which remains well above the 
room temperature Boltzmann limit of 60 mV/decade. 

A schematic of the BLG device configuration is shown in Fig. [ija) where a BLG is sand- 
wiched between a top and a bottom gate stack. Each gate stack consists of a gate dielectric 
with an equivalent oxide thickness (EOT) of lnm and a gate contact metal. We assume per- 
fect dielectric-BLG interface and zero flatband voltage. The voltage at the top gate (Vt) and 
the bottom gate (V&) can be varied independently. The self-consistent electronic structure of 
this device configuration is determined by the bandstructure of the BLG using Tight Bind- 
ing Method ([HO [TO]) coupled with the 1-D Poisson equation. Taking the center of the two 
graphene layers at z = 0, the charge density is given by p(z) = qn{z) where q is the electronic 
charge and n(z) is obtained as a subtraction between the hole (n^(z)) and electron (n e (z)) 
carrier density: 
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where f(Eij(k)) is the Fermi-Dirac probability of the state ({i,j},k) at temperature T. i and 
j are the valence and conduction band indices respectively. The energy eigenvalues Eij(k) are 
obtained from the tight binding bandstructure taking only p z orbital into account, with an 
intra-layer overlap integral (5) of 0.129, the intra-layer hopping t as — 3.033eV and inter-layer 
hopping t± as — 0.365eV pQ. To obtain the wavefunction ip!?(z), we assume normalized Gaussian 
orbital as the basis function. The wavefunctions are set to zero at the BLG-dielectric interfaces. 

We validate the above self-consistent method by comparing the bias dependent bandgaps 
of BLG experimentally obtained in [6]. In this case, Vt is set to — Vb which breaks the inversion 
symmetry of the BLG opening a bandgap [HE]. The results are shown in Fig. [TJb) where the 
bandgap (E g ) is plotted as a function of the average displacement vector D, defined in the same 
way as in [6]. Here the relative permittivity (e r ) of BLG is varied as a parameter to find that 
e r =1.8 provides the best fit to the experimental data and the same value is used in the Poisson 
equation to generate the results discussed in the rest of the paper. Fig. []Jc) shows the carrier 
density profiles along z at Z}=3V/nm. This clearly shows a strong charge polarization in the 
bilayer, though the system as a whole remains almost intrinsic for Vt = —Vb- 

The recent development of significant tunability of bandgap in a bilayer graphene ( [5]- 110| ) 
brings the possibility of significant reduction of off current in a BLG transistor. In the following, 
we show that the bias dependent electronic structure and the corresponding self-consistent 
electrostatics play a major role in determining the intrinsic limits of the subthreshold slopes 
and the on/off current ratio in an ideal long channel BLG MOSFET. 

Before discussing the self-consistent results, we first show that even in absence of the screen- 
ing effect, the intrinsic bias dependence of the electronic structure in BLG can be a major cause 
for subthreshold slope degradation. This arises from the dependence of the rate of change of 
electron and hole barrier height (B e h) with gate bias. This rate, ideally, should be as high as 
possible, with the maximum possible value being 1, which corresponds to a subthreshold slope 
of 60mV/decade. This is obtained by finding the roots A of the secular equation of the bilayer 
Hamiltonian matrix (with 5 = 0): 

[04-A) 2 -t 2 | Sfc | 2 ] l( Vb -\) 2 -t 2 \s k \ 2 ] 

+tl(Vt-X)(V b -X) = (2) 
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where Sk = 1 + e^-" 1 + e 1 ^' 0-2 , d\ and d% are lattice vectors in graphene. In Fig. [l|d), we 
plot the rate of change of B e for different directions in the (Vt, Vb) space in polar coordinates 
with 9 = tan~ 1 (Vt/Vb). This clearly shows that the rate is maximum (=1) along Vt = Vb + Vo, 
which is essentially the constant bandgap locus. As we deviate from this direction, the rate 
degrades from 1 with minimum along Vt = —Vb, leading to the poor subthreshold slope along 
anti-symmetric bias case. As observed from Fig. [l^d), this is a fairly generic result for a large 
variety of points in the subthreshold region. 

We now consider a chemically undoped long channel BLG MOSFET with metal source and 
drain having low drain bias. The drain current is assumed to be completely dominated by 
the thermionic carriers, neglecting tunneling. Under these approximations, the electron barrier 
height B e is plotted in Fig. [2^a) as a function of Vt and Vb- In the same plot, we also show a 
number of possible paths p from transistor off state (labeled X p points) to on state (labeled Y p 
points). 

Due to the ambipolar nature of the device, it is important to select the appropriate off state 
for the device. The total integrated carrier density (N = N e + Nh) is shown in Fig. [2](b) in 
the (Vt,Vb) space to help choose the appropriate set of off state points. In Fig. [2|a), Xj.,4 
corresponds to the point of maximum bandgap in that range. X2 corresponds to the maximum 
electron barrier height though it reduces the hole barrier height, hence increasing N. X3 is 
an intermediate point with symmetric electron and hole barrier heights. We now compare the 
characteristics among these off to on paths shown in Fig. [2^a), namely (1): Vb is fixed at IV 
and Vt is varied from —IV to IV, (2): Vb is fixed at 0.5V and Vt is varied from —IV to IV, 
(3): Vb = Vt + IV, and Vt is varied from —IV to 0V, which is an almost constant bandgap 
operation (with E g ~ 0.17eV), (4): Vb = —Vt and Vt is varied from —IV to 0V, which is the 
maximum bandgap tunability path. It should be noted that (1) and (2) correspond to single 
gate operation with back gate at fixed bias. On the other hand, (3) and (4) correspond to double 
gate operation, although unconventional, with top and bottom gates at different biases. The 
conventional double gate operation with Vt = Vf, is of little significance due to zero bandgap. 
Fig. [2^c) shows N as a function of Vt in these four cases. It is clearly observed that the best 
on/off ratio obtained is ~100. Also, path (1) and (4) will have lower off current than (2) and 
(3), which is also expected from the Fig. pfb). Path (4) is found to have low on/off ratio 
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due to the poor tunability of barrier heights (though bandgap tunability is maximum) in the 
anti-symmetric case. 

To get insights into the above, we find the drain current in a long channel transistor as 
Id = qNr, where r is the transport factor. In a MOSFET configuration, r is a weak function 
of gate voltage. Hence, the subthreshold slope is SS = ln(10) gj^} d ) ~ £ra(10) g^m) ■ ^ n Fig- 
[2](d), the computed subthreshold slopes are plotted for the paths (1) to (4) with varying Vt. 
Note that, unlike conventional MOSFET, the computed subthreshold slope is not independent 
of gate bias in the BLG transistor. 

Analytically, we can convert the summation used in Eq. [T] into an integral to obtain N: 

(3) 



N 



D v (E)(l - f{E))dE + / D c {E)f(E)dE 
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where only one band is assumed to contribute in the conduction band (CB) and the valence 
band (VB). E v>c are the VB max and CB m j„ obtained from the self-consistent electrostatics. 
D v c are the 2-D DOS of VB and CB respectively, and we replace them by an average DOS D 
using the fact that close to band edge, D V)C change much slowly as compared to Fermi function 
[9]. Utilizing CB-VB symmetry, we find, 



N = 2k B TD\n 



(l + e-^)(l + e-^) (4) 



where B' e = ^ = =fg^ and B h = ^ = ^=^. As D is weakly dependent on V t [9j, we 
obtain 

SS = -K10), feTKMM (6) 

where k = (1 + e~ B «)(l + e~ B 'h) and E' g = ^ = Be k ^ h . In a conventional Si MOSFET, in 
the subthreshold region, B' e , B' h » 0, E g is constant and no gate field screening occurs due to 
small carrier density (I q Qy^ l = l)- Under these conditions, it is straight forward to show that 
SS reduces to ln(10)(-|^) (=60mV/decade at T = 300K). Unfortunately, for bilayer graphene, 
these three conditions are never satisfied simultaneously which degrades the subthreshold slope. 
To open up a bandgap in BLG, unlike conventional semiconductor, we need to apply significantly 
large fields at the two gates. This in turn causes a large number of carriers to be present in the 
device, even in subthreshold (though it is possible that the overall device is close to neutral, 
like in the case of Vt = — Vj,, Fig. [jjc)). These carriers cause a strong screening of the gate field 



and hence, a reduction in the derivatives in the denominator of Eq. [5] Thus, even if we create 
a significant bandgap using vertical field, the carriers thus generated degrade the subthreshold 
slope. 

In the anti-symmetric bias case (path (4)), using B e = Bh = we find from Eq. [5] that the 
SS limit reduces to 60/|^||. As |Jgpf | is intrinsically small in this direction (see Fig. |T|(d) ) , 
SS degrades significantly. On the other hand, along the constant bandgap paths (e.g. path 
3), the SS limit again reduces to 60/|4jp||. However, in this case, the strong screening due to 
large number of carriers reduces the denominator, limiting SS much above Boltzmann limit. 
Finally, we have numerically performed a global search in the (Vt,V&) space with \Vtb\ < IV, 
which shows that the best SS for all possible constant bandgap paths is around 200mV/decade, 
whereas along all possible constant V& paths is around 275mV/decade. 

Thus, we find that along the constant bandgap loci, the fundamental limit of SS arising 
from the bias dependent electronic structure is 60mV/decade (see Fig. [ijd)), but the screening 
due to the presence of large number of carriers arising from the unavoidable gate field degrades 
the subthreshold slope significantly. However, as the operating condition moves away from 
this direction, both the electronic structure as well as the screening effect play major roles 
to degrade SS. Nonetheless, the shown subthreshold slope limits are significantly less than 
the most recent experimental data reported [11] where the extracted value of the minimum 
subthreshold slope is ~550 mV/decade. This clearly shows that there is room for significant 
technological improvement in bilayer graphene, including efforts in improving transport derived 
bandgap, effective oxide thickness and choosing the optimized path from on state to off state. 

In conclusion, in this work, we have developed a coupled self-consistent Poisson-bandstructure 
solver for dual gated bilayer graphene. The bias dependent bandgap results obtained from the 
solver have been benchmarked against published experimental data. Finally, it has been shown 
that the bias dependent electronic structure and the self-consistent electrostatics limit the sub- 
threshold slope obtained in such a transistor well above the Boltzmann limit of 60mV/decade, 
but again much below the results experimentally shown till date. 
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Figure 1: (a): Schematic of a bilayer graphene (BLG) sandwiched between two gate stacks, 
(b): Verification of numerically obtained bandgap results with experimental data in [6], shown 
as red dots, (c): Carrier density profiles from top layer to bottom layer, (d): Polar plot 
of the rate of change of B e (normalized by q) along different directions using Eq. [2j (with 
angle 6 = ian _1 (Vt/Vb)) in the (Vt, V&) space for two sets of subthreshold points: i) Vt=— IV, 
V b £ [0V,1V] and ii) V 4 =-0.5V, V b G [OV, 0.5V]. All the curves show similar directional 
dependence. 
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Figure 2: (a): Electron Barrier height B e as a function of Vt and Vb- The four different off 
to on state arrow paths are numbered as (1) to (4). (b): Integrated carrier density (N) as a 
function of Vt and (c): N along the different off to on paths (d): The absolute values of the 
computed subthreshold slopes for the paths (l)-(4) without any non-ideality. 
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